The 22nd Texas Symposium on Relativistic Astrophysics, Stanford University, December 2004 



1 



Gravitational Radiation from Coalescing Supermassive Black Hole 
Binaries in a Hierarchical Galaxy Formation Model 

Motohiro Enoki * \ Kaiki T. Inoue 2 , Masahiro Nagashima 3 and Naoshi Sugiyama 1 

1 National Astronomical Observatory of Japan, Tokyo 181-8588, Japan 

2 Department of Natural Science and Engineering, Kinki University, Higashi-Osaka, Osaka, 577-8502, 
Japan and 

3 Department of Physics, Kyoto University, Sakyo-ku, Kyoto, 606-8502, Japan 



in 
o 
o 

(N 

X> 

<D 

oo 

(N 
> 

o\ : 

(N . 

in ■ 

(N ■ 

o : 

in . 
o ■ 

Or 

6 ■ 
a: 



We investigate the expected gravitational wave emission from coalescing supermassive black hole (SMBH) 
binaries resulting from mergers of their host galaxies. We employ a semi-analytic model of galaxy and quasar 
formation based on the hierarchical clustering scenario to estimate the amplitude of the expected stochastic 
gravitational wave background owing to inspiraling SMBH binaries and bursts rates owing to the SMBH binary 
coalescence events. We find that the characteristic strain amplitude of the background radiation is h c [f) ~ 
10- ie (//lAtHz)' 2 / 3 for / < 1/xHz. The main contribution to the total strain amplitude of the background 
radiation comes from SMBH coalescence events at < z < 1. We also find that a future space-based gravitational 
wave interferometer such as the planned Laser Interferometer Space Antenna (LISA) might detect intense 
gravitational wave bursts associated with coalescence of SMBH binaries with total mass 

Aftot < 1O 7 M at 

z 2 at a rate ~ 1.0 yr" 1 . Our model predicts that burst signals with a larger amplitude ftburst ~ 10~ 15 
correspond to coalescence events of massive SMBH binary with total mass Aftot ~ 10 8 Mq at low redshift z ^ 1 



at a rate ~ 0.1 yr 1 whereas those with a smaller amplitude h^ ur 



10" 



of less massive SMBH binary with total mass Aft, 



10 6 M Q at high redshift z > 3. 



correspond to coalescence events 



1. Introduction 



In recent years, there has been increasing observa- 
tional evidence that many nearby galaxies have cen- 
tral supermassive black holes (SMBHs) in the mass 
range of 10 6 — 10 9 Af Q , and that their physical prop- 
erties correlate with those of spheroids of their host 
galaxies. This suggests that the formation of SMBHs 
physically links to the formation of spheroids that har- 
bor the SMBHs. Therefore, in order to study the for- 
mation and evolution of SMBHs, it is necessary to 
construct a model that includes galaxy formation pro- 
cesses. 

In the standard hierarchical structure formation 
scenario in a cold dark matter (CDM) universe, dark- 
matter halos {dark halos) cluster gravitationally and 
merge together. In each of merged dark halos, a 
galaxy is formed as a result of radiative gas cool- 
ing, star formation, and supernova feedback. Several 
galaxies in a common dark halo sometimes merge to- 
gether and a more massive galaxy is assembled. When 
galaxies merge, SMBHs in the centers of the galaxies 
sink toward the center of the new merged galaxy and 
form a SMBH binary subsequently. If the binary loses 
enough energy and angular momentum, it will evolve 
to the gravitational wave emitting regime and begin 
inspiraling, eventually coalesces with a gravitational 
wave burst. 

An ensemble of gravitational waves from a number 
of inspiraling SMBH binaries at different redshift can 



be observed as a stochastic background at frequencies 
~ In — 1/iHz, which can be detected by pulsar tim- 
ing measurements. Future space interferometers such 
as the Laser Interferometer Space Antenna (LISA) 
might detect quasi-monochromatic wave sources as- 
sociated with inspiraling SMBH binaries and gravi- 
tational wave bursts associated with SMBH binary 
coalescence [e.g. Q. 

To date, a number of attempts have been made 
to calculate the SMBH coalescence rate. Some au- 
thors use phenomenological models of galaxy merg- 
ers based on number counts of quasars and spheroids 
[e.g. 0, 0. Others use merger rates of dark halos 
(not galaxies) [e.g. 01- However, none of these models 
include baryonic gas evolution and galaxy formation 
processes. Because SMBH formation process is rele- 
vant to spheroids of host galaxies rather than to dark 
halos, we need to evaluate how the baryonic gas pro- 
cesses such as star formation and supernova feedback 
affect the SMBH formation process. 

In this study, we estimate the SMBH coalescence 
rate using a new semi-analytic (SA) model |4| (an ex- 
tended model 5]) in which the SMBH formation is 
incorporated into the galaxy formation. Then, we cal- 
culate the spectrum of gravitational wave background 
from inspiraling SMBH binaries, based on the formu- 
lation given by Jaffe & Backer |2| and we compare our 
result with that from a pulsar timing measurement. 
We also estimate the event rate of gravitational wave 
bursts from SMBH coalescence events that might be 
detected by future planned space laser interferome- 
ters, based on an argument in The details of for- 
mulation, results, discussion and references are given 
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2. Galaxy Merger / Black Hole 
Coalescence Rate 

Here we briefly describe our SA model for galaxy 
formation and SMBH growth. The details are shown 
in @ and 0. 

2.1. The model of galaxy formation 

First, we construct Monte Carlo realizations of 
merging histories of dark halos from the present to 
higher rcdshifts. Merging histories of dark halos de- 
pend on the cosmological model. The adopted cos- 
mological model is a low-density, spatially flat cold 
dark matter (ACDM) universe with the present den- 
sity parameter, Q m = 0.3, the cosmological con- 
stant, S1a = 0.7, the Hubble constant h = 0.7 
(h = i?o/100 km s^ 1 Mpc" 1 ) and the present rms 
density fluctuation in spheres of 8/i _1 Mpc radius, 
cj 8 = 0.9. The highest redshift in each merging path 
which depends on the present dark halo mass, is about 
z ~ 20 - 30. 

Next, in each merging path of dark halos, we cal- 
culate the evolution of the baryonic component from 
higher redshifts to the present using simple analytic 
models for gas cooling, star formation, supernova feed- 
back, galaxy merging and other processes. When a 
dark halo collapses, the gas in the halo is shock-heated 
to the virial temperature of the halo (the hot gas). At 
the same time, the gas in dense regions of the halo 
cools owing to efficient radiative cooling and sinks to 
the center of the halo and settle into a rotationally 
supported disk until the subsequent collapse of the 
dark halo. We call this cooled gas the cold gas. Stars 
are formed from the cold gas. With star formation, 
supernovae occur and heat up the surrounding cold 
gas to the hot gas phase (supernova feedback). 

When several progenitor halos have merged, a 
newly formed larger dark halo contains at least two 
or more galaxies which had originally resided in the 
individual progenitor halos. We identify the cen- 
tral galaxy in the new common halo with the cen- 
tral galaxy contained in the most massive of the pro- 
genitor halos. Other galaxies are regarded as satel- 
lite galaxies. After each merging of halo, these satel- 
lites merge by either dynamical friction or random 
collision. Satellite galaxies merge with the central 
galaxy in the dynamical friction timescale Q. Satel- 
lite galaxies sometimes merge with each other in the 
timescale of random collisions. Under the condition 
that the satellite galaxies gravitationally bound and 
merge during encounters, this timescale is derived by 
Makino & Hut |g. If the mass ratio, / = mijmi, 
is larger than a certain critical value of /bulge, we as- 
sume that a starburst occurs, and that all of the cold 
gas turns into stars and hot gas, which fills the dark 
halo, and all of the stars populate the bulge of a new 



galaxy. On the other hand, if / < /bulge, no starburst 
occurs and a smaller galaxy is simply absorbed into 
the disk of a larger galaxy. 

Model parameters are determined by a comparison 
with observations of galaxies in the local Universe, 
such as luminosity functions and the cold gas mass 
fraction in spiral galaxies. Our SA model can repro- 
duce galaxy number counts and photometric redshift 
distribution of galaxies in the Hubble Deep Field 

2.2. The growth of SMBHs 

In our model, it is assumed that SMBHs grow by 
coalescence when their host galaxies merge and are 
fueled by accreted cold gas during major mergers of 
galaxies. When the host galaxies merge, pre-existing 
SMBHs sink to the center of the new merged galaxy 
due to dynamical friction (or other mechanisms such 
as gas dynamics), evolve to the gravitational wave 
emission regime and eventually coalesce. Although 
the timescale for this process is unknown, for the 
sake of simplicity we assume that SMBHs instan- 
taneously evolve to the gravitational wave emission 
regime and coalesce. Gas-dynamical simulations have 
demonstrated that the major merger of galaxies can 
drive substantial gaseous inflows and trigger starburst 
activity. Thus, it is reasonable to assume that dur- 
ing a major merger of galaxies, a certain fraction of 
the cold gas that is proportional to the total mass 
of stars newly formed at starburst accretes onto the 
newly formed SMBH and this accretion process leads 
to a quasar activity. Under this assumption, the mass 
of cold gas accreted on a SMBH is given by 

-Mace = /BHAAf^burst, (1) 

where /bh is a constant and AAf*. burst is the total 
mass of stars formed at starburst. The free param- 
eter of /bh is fixed by matching the observed rela- 
tion between a spheroid mass and a black hole mass 
M B H/M sp heroid = 0.001-0.006 (e.g. @); we find that 
the favorable value of /bh is nearly 0.03. 

Figure ^( a ) shows the black hole mass functions in 
our model at a series of various redshifts. In this fig- 
ure, we superpose the observed black hole mass func- 
tion at z = [3. The predicted mass function is 
quite consistent with the observation. Our galaxy for- 
mation model includes dynamical friction and random 
collision as galaxy merging mechanisms. For compar- 
ison, in figure H(b), we also plot the black hole mass 
functions at z = of other two models: no random 
collision model (no rc model) and no dynamical fric- 
tion model (no df model) . In the no rc model and the 
no df model, mergers owing to random collision and 
dynamical friction, respectively, are switched off. This 
figure shows that the mass function for low mass black 
holes are determined by random collisions between 
satellite galaxies and that for high mass black holes 
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Figure 1: (a) Black hole mass function of the model at a series of redshifts. The symbols with errorbars are the 
observed black hole mass function at z = E3. (b) Black hole mass function at z = for different galaxy merger 
models, which are both model, no random collision model (no rc model) and no dynamical friction model (no df 
model). Both model includes dynamical friction and random collision as the galaxy merging mechanism. 



are influenced by dynamical friction. The shape of 
black hole mass function depends on detailed gas pro- 
cesses. The important contribution of mass increment 
of SMBHs in central galaxies is the cold gas, which 
is accreted to only central galaxies. SNe feedback re- 
moves this cold gas more efficiently in smaller galaxies 
with Vdrc < 280km s _1 (M gal < 10 12 M@). Thus, the 
growth of the SMBHs in small galaxies suffers from 
SNe feedback. In the no rc model, SMBHs mainly ex- 
ist central galaxies. Therefore, the shape of the black 
hole mass function in the no rc model has a bump at 
high mass end (Mbb ~ 1O 9 M ) which corresponds to 
SMBHs in the central galaxies (M gal ~ 10 12 Af o ). On 
the other hand, in the no df model, high mass SMBHs 
cannot be produced since galaxies cannot merge with 
the massive central galaxy. 

Using the SA model incorporated with this SMBH 
growth model, we estimate the comoving number 
density, n c (M\,M 2 , z)dM\dM 2 dz, of the coalescing 
SMBH binaries with mass M x ~ M\ + dM\ and 
M 2 ~ M 2 + dM 2 at z in a given redshift interval dz. 



3. Gravitational Radiation 

3.1. Background radiation from SMBH 
binaries 

In order to calculate the spectrum of gravitational 
wave background radiation from SMBH binaries, we 
improve the formulation of Jaffe & Backer 0|. The 
details are described in |fj . The spectrum of the grav- 
itational wave background radiation which we finally 
obtain is 

hiu) = Jdz dMl d M 2 (^|^) 5/3 wr 4/3 

x(l + «)- 1 /3„ c (M 1 , M 2 , z) 0(/ max - /), (2) 



where M chirp = [M 1 M 2 (M 1 + Af 2 )~ 1/3 ] 3/5 is the chirp 
mass of the system and c is the speed of light, / 
is the observed frequency of the gravitational wave 
from the binary in a circular orbit, the / max is max- 
imum frequency and 0(x) is the step function. As a 
SMBH binary evolves with time, the frequency be- 
comes higher. We assume that the binary orbit is 
quasi-stationary until the radius equals to 3i?s, where 
Rs is the Schwarzschild radius : the radius of the in- 
nermost stable circular orbit (ISCO) for a particle and 
a non-rotating black hole. Then / max is 

c 3 / Mo \ 1/2 

/,„«(M„M 2 ,„ = p75 _^- 5 (l + _) 

/ TUT \ 1 / 2 

x(l + B ?) Hz, (3, 

where M x and M 2 are SMBH masses (Mi > M 2 ). 

As shown in figure[21 the spectrum changes its slope 
at / ~ 1/iHz owing to lack of power associated with 
the upper limit frequency, / max . This feature is con- 
sistent with the results of Sesana et al. 01 • The pre- 
dicted strain spectrum is h c (f) ~ 10~ 16 (//l/iHz)~ 2 / 3 
for / < 1/iHz, just below the current limit from the 
pulsar timing measurements |T^ . In our model, we 
assume that SMBHs coalesce simultaneously when 
their host galaxies merge. Therefore, the efficiency 
of SMBH coalescence is maximum and the predicted 
amplitude of gravitational wave spectrum should be 
interpreted as the upper limit. 

In figure |21 (a), we plot the spectra from binaries in 
different redshift intervals. This figure shows that the 
total spectrum of background radiation comes from 
coalescing SMBH binaries at low redshift, < z < 1. 
Although the coalescence rate at low redshift is lower 
than at high redshift, the main contribution to the 
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Figure 2: Spectrum of gravitational wave background 
radiation, h c (f), from SMBH binaries in different 
redshift intervals (a) and in different total mass intervals 
(b). The filled reverse triangle in (a) shows the current 
limit from pulsar timing measurements |T^|. 



background radiation is the coalescing SMBH bina- 
ries at low redshift, < z < 1. This is because the 
distance from an observer to the SMBH binaries at 
low redshift is shorter and the mass of SMBHs at 
low redshift is higher. In figure l^b), we also plot 
the spectra from binaries in different total mass inter- 
vals (Mtot = Mi + M2). This figure shows that for 
/ > 10~ 4 Hz the total spectrum of background radia- 
tion comes from coalescing SMBH binaries with total 
mass M tot ^ 10 8 M Q . 

When SMBHs are spinning and/or when their 
masses are comparable, the definition of ISCO be- 
comes vague and our assumption that the cutoff fre- 
quency / max corresponds to 3i?s may not be correct. 
To see the effects of the cutoff frequency, we plot the 
spectra for different values of / max , corresponding to 
3i?s, 30i?Sj and no frequency cut off, respectively (fig. 
(a)). Lowering / max causes a suppression in the 
stochastic background at high frequencies / < 10 _7 Hz 
since a large portion of high frequency modes are cut 
off. 

Our galaxy formation model incorporates dynam- 
ical friction and random collision as galaxy merging 
mechanisms. In order to examine the effect of these 
two galaxy merger mechanisms on the spectrum of 
gravitational wave background radiation, in figure 
(b), we also plot the spectrum of background radia- 
tion of other two models: no rc model and no df model. 
The no df model can not produce higher mass SMBH 
(see fig. Bib)). Consequently, the spectrum in no df 
model bends at larger frequency since the number of 
SMBHs with smaller / max decreases. Furthermore, in 
the no df model, the amplitude of gravitational waves 
from each binary becomes smaller . However, the coa- 
lescence rate in no df model is higher than the rate in 
no rc model. As a result, the amplitude of the spec- 
trum in no df model is roughly equal to the amplitude 
in no rc model. 



Figure 3: Spectrum of gravitational wave background 
radiation, h c (f), in different / max (a) and in different 
galaxy merger models (b). 



3.2. Gravitational wave bursts from 
SMBH coalescence 

After an inspiraling phase, SMBHs plunge into a 
final death inspiral and merge to form a single black 
hole. We call a set of a plunge and a subsequent early 
non-linear ring-down phase as a burst. We estimate 
the expected burst event rate per observers' time using 
the amplitude of burst gravitational wave given in U 
and the SMBH coalescence rate calculated by our SA 
model, n c (Mi, M 2 , z). 

The characteristic amplitude of the burst gravita- 
tional wave which we finally obtain is 



3 3/4 e l/2 GM 



"burst 



2 1 / 2 irc 2 D(z) 



7.8 x 1(T 16 



Get) 



1/2 



M t „ 



10 s Mq 



D{z) 



1 Gpc 



(4) 



where e is the efficiency of the energy release, D(z) is 
the comoving distance to the source from the observer 
and f c is the observed characteristic frequency. f c 
from the gravitational wave burst occurring at z is 



fc = 



3 3 / 2 GM to t (1 + z) 
M tot 



3.9 x 1(T 4 



10 8 Mr. 



(l + z)- 1 Hz. (5) 



Then, we obtain n hmst (h hmst , f c , z) dh hul . st df c dz , 
which is the comoving number density of gravitational 
wave burst events occurring at z in a given redshift in- 
terval dz with amplitude /iburst ~ Tiburst + rf^burst and 
with characteristic frequency f c ~ f c + df c . Therefore, 
the expected event rates of gravitational wave bursts 
per observers' time with amplitude /iburst ~ ^burst + 
dhburat and characteristic frequency f c ~ f c + df c , 
z-Wst (iburst, fc) rf/iburst df c is given by 



^burst (^burst , fc 



"burst (iburst, fc, z)^-dz,(6) 

dt 
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Figure 4: Integrated event rate of gravitational wave bursts from SMBH coalescence per observers' time unit a year, 
!-'burst[log(/iburst)], in different redshifts (a) and in different total mass ranges (b). 




log(f c [Hz]) log(fc [Hz]) 

Figure 5: Expected signals of gravitational wave bursts from SMBH coalescence, (a) The horizontally hatched area 
shows the region, ^burst[log(/iburst), log(/ c )] > lyr -1 . The solid curve indicates the instrumental noise threshold for one 
year of LISA observations, (b) The vertically hatched area shows the region, fburst > l/5yr _1 and the diagnal 
cross-hatched area show the region, fburst > 3yr _1 . The dot-dashed and the short dashed lines indicates the 
instrumental noise threshold for 5 year and 1/3 year of LISA observations, respectively. 



where to is the observers' time and dV is the comov- 
ing volume element at z. The integrated event rates of 
gravitational wave bursts per observers' time with am- 
plitude /iburst ~ ^burst +<^burst) ^burst (/iburst) d/lburst i 

are given respectively as follows, 

Z4)urst (/iburst ) = / ^burst (/iburst, fc)df c . (7) 

In figure we plot the total integrated event rates 
of gravitational wave bursts and integrated event rates 
in different redshift intervals (fig. EI a )) and in dif- 
ferent total mass intervals in (fig. Bib)). Here we 
set the efficiency of the energy release e = 0.1, while 
the precise value of this parameter is unknown. In 
most typical events, a conversion efficiency will prob- 
ably be a few percent. From equation (QJ, one can 
see that the change of efficiency results in the par- 
allel displacement in the horizontal direction in fig- 
ure 0] The shape of i^w-st (/iburst) reflects the black 
hole mass functions and the SMBH coalescence rates, 
which depend on the complex galaxy formation pro- 



cesses. In figure 0] (a), one can notice that there are 
two peaks in the event rate in terms of /iburst- A peak 
at /iburst ~ 1CP 17 corresponds to bursts from SMBH 
binaries with M to t < 1O 6 M0 whose total number is 
the largest at high redshift z > 3. Another peak at 
/iburst ~ 10~ 15 corresponds to bursts from SMBH bi- 
naries with 10 7 M Q < M to t < 10 8 M Q whose coales- 
cence probability is the largest at low redshift z < 3. 
Figure 0] (b) indicates that burst signals with large 
amplitude (/iburst ^ 10~ 15 ) correspond to coalescence 
of "massive" SMBH binaries with M tot > 10 8 M Q oc- 
curring at z ^ 1. This is because the distance from 
the earth to SMBHs at low redshift is shorter and 
the mass of SMBHs at low redshift is larger. On 
the other hand, burst signals with small amplitude 
(/iburst ^ 10~ 17 ) corresponds to coalescence events of 
"less massive" SMBH binaries with M tot < 10 7 M© 
occurring at z > 2. These events dominate the ex- 
pected burst event rate provided that the sensitivity 
of the detector is sufficiently good. This feature is 
quite important because it breaks the degeneracy be- 
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tween mass and distance. 

Figure [S] shows that the expected region for sig- 
nal of gravitational wave bursts and the instrumen- 
tal noise threshold for LISA, hi nst . We compute /linst 
from the fitting formula for the spectral instrumental 
noise density of LISA 13] . The expected region for 
^burst [log(/iburst),log(/ c )] > 1 yr 1 is above this in- 
strumental noise threshold. For comparison, we show 
the region for ^burst[fog(ftburst), fog(/c)] > 1/5 yr" 1 
and «4>urst[log(/iburst))log(/ c )] > 3 yr" 1 in figureEl(b). 

From figure 0] and 03 we conclude that the LISA can 
detect intense bursts of gravitational waves at a rate of 
~ l.Oyr" 1 assuming that dominant part of these burst 
events occur at z > 2. Even in the case of e = 0.001, 
the LISA can detect intense bursts of gravitational 
waves in one year observation, since ft-burst oc e 1 / 2 . In 
addition, we find that large amplitude /iburst ~ 10 -15 
signals correspond to coalescence events of massive 
SMBH binaries M tot ~ 1O 8 M at low redshift z < 1 
and small amplitude /iburst ~ 10~ 17 signals corre- 
spond to less massive SMBH binaries M to t ~ 10 6 Mq 
at high redshift z > 3. 

4. Summary & Conclusions 

In this study, we have estimated the coalescence 
rate of SMBH binaries in the centers of galaxies using 
a new S A model of galaxy and quasar formation £| . 
Then, we calculated the spectrum of the gravitational 
wave background from inspiraling SMBH binaries and 
estimated the expected amplitudes and event rates of 
intense bursts ofgravitational waves from coalescing 
SMBH binaries @. 

Our SA model includes dynamical friction and ran- 
dom collision as galaxy merging mechanisms, and as- 
sumes that a SMBH is fueled by accretion of cold 
gas during a major merger of galaxies leading to a 
spheroid formation, and that SMBHs coalesce simul- 
taneously when host galaxies merge. Many previ- 
ous other studies have paid attention to only SMBH 
growth and did not take galaxy formation processes 
into account. For investigating the relations between 
SMBH growth and galaxy formation processes, SA 
methods of gala xy and SMBH formation are suitable 
tools [e.g. Hill . 

We have found that the gravitational wave back- 
ground radiation spectrum for / < 1/iHz has a char- 
acteristic strain h c (f) ~ 10" 16 (//l^Hz)" 2 / 3 just be- 
low the detection limit from the current measurements 
of the pulsar timing. The slope of the spectrum for 
/ ^ 1/xHz gets steep owing to the upper limit in fre- 
quency set by the radius of the ISCO. The stochastic 
background radiation mainly comes from inspiraling 
SMBH binaries at < z < 1. Therefore, the back- 
ground radiation can probe inspiraling SMBH binaries 
at low redshift . 

We have also found that LISA might detect intense 



bursts of gravitational waves owing to the SMBH co- 
alescence events at a rate 0.1 ~ 1.0 yr -1 and that the 
main contribution to the event rate comes from SMBH 
binary coalescence events at high redshift z > 2. Our 
model predicts that burst signals with a large ampli- 
tude correspond to coalescence of large mass SMBH 
binaries at low redshift while those with a small ampli- 
tude correspond to coalescence of small mass SMBH 
binaries at high redshift. This prediction can be tested 
by future measurements of the amplitude and the 
phase evolution in gravitational waves from inspiral- 
ing SMBH binaries [13j . Comparing these predictions 
with observations in future, we can put a stringent 
constraint on SMBH formation and evolution models. 
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